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ABSTRACT 

We study Type I migration of a planet in a radiatively efficient disk using global 
two dimensional hydrodynamic simulations. The large positive corotation torque is 
exerted on a planet by an adiabatic disk at early times when the disk has the steep 
negative entropy gradient. The gas on the horseshoe orbit of the planet is compressed 
adiabatically during the change of the orbit from the slow orbit to the fast orbit, 
increasing its density and exerting the positive torque on the planet. The planet would 
migrate outward in the adiabatic disk before saturation sets in. We further study 
the effect of energy dissipation by radiation on Type I migration of the planet. The 
corotation torque decreases when the energy dissipates effectively because the density 
of the gas on the horseshoe orbit docs not increase by the compression compared 
with the gas of the adiabatic disk. The total torque is mainly determined by the 
negative Lindblad torque and becomes negative. The planet migrates inward toward 
the central star in the radiatively efficient disk. The migration velocity is dependent on 
the radiative efficiency and greatly reduced if the radiative cooling works inefficiently. 

Key words: hydrodynamics- radiative transfer-method: numerical-type I migration: 
protoplanet- planetary system: gravitational interaction 



1 INTRODUCTION 

The density waves in a protoplanetary disk are ex- 
cited by the gravitational inter actions with a planet 
jGoldreich fc Tremainelll979l . ll980T ). The tidal torque acts 
from the disk on a planet and generally leads to the migra- 
tion of the planet. This process is known as Type I migration 
and plays an important role in the formation of planets. It is 
suggested that Type I migration is crucial to und erstand the 
semim ajor axis distribution of extrasolar planets (|lda fc LirJ 
|200S| ). 

There are two kinds of torques acting on a planet: one 
is the Lindblad torque due to spiral arms in a disk and the 
other is the corotation torque due to the gas which peri- 
odically exchanges angular momentum with a planet on its 
horseshoe orbit. The Lindblad torque occurs at the Lindblad 
resonances a nd leads to the inward migration of a planet 
l|Wardl fl997l ). On the other hand, the corot ation torque 
is dependent on the physical state of a disk (|Wardl [l99ll: 
iPaardekooper fc Papaloizoull2009bl ) and invokes the inward 
or outward migration of a planet. The sum of the two torque 
values determines the direction and magnitude of migration 
of a planet. 

Type I migration was studied analytically and numeri- 
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cally using isothermal dWardlll986l Il989l: [Artvmowicj[l993 : 
iTanaka et all 120021: ICasoli fc Massed 120091') and adiabatic 



disks ( Paai dek ooper fc Mc llcma 2006b: iBaruteau fc Massct 
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20091 : IPaardekooper et al.ll2010T ) . It was shown that a planet 
in an isothermal disk migrates inward toward a central star 
in a shorter time than a lifetime of a disk jTanaka et al.l 
120021 ). The rapid inward migration suggests that planets in 
an isothermal disk cannot survive or stay in a far region of a 
disk. In an adiabatic disk with the negative gradient of the 
entropy, the density of the gas on the horseshoe orbit of a 
planet increases when the gas approaches and moves away 
from the planet, yielding the positive corotation torque. The 
positive corotation torque might become large enough to 
chang e the direction of the migr ation of a planet to be out- 
ward (|Baruteau fc Massetl 12001 ). 

Isothermal and adiabatic disks are idealized disks and 
different from real disks with dissipative processes. When 
a planet excites density waves in a disk by the grav- 
itational interactions, the density waves can be altered 
by dissipative processes. Numerical simulations and the 
linear calculation were carried out, taking into account 
viscosity and thermal conduction as dissipative processes 



Yu et all 



in disks ( Paard ekooper 



201011. 



lduction as dissipati ve processes 
Papaloizoul 12008 : iLi et all 120091 : 



With the use of the linear calculation, 
Muto fc Inutsukal (|2009l ) examined the one side torque act- 
ing on a planet, including viscosity expressed by the a 
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prescription. It was found that the one side torque does 
not change greatly unless the value of a is larger than 
0.01, which happens to be the upper limit of a obtained 
from observations of disks l|Hartmann et al1ll998l ). This is 
not important in laminar disks and is not included in this 
study.However, viscosity does play a role in sustaining the 
corotation torque over longer time scales, which we do not 
consider in this paper. IPaardekooper fc Papaloizoul (|2008T ) 
included thermal diffusion in a disk as another dissipative 
process. They showed that the corotation torque decreases 
with an increase in the thermal diffusion and the total torque 
reduces to that in an isothermal disk. A planet migrates in a 
different direction, depending on the strength of the thermal 
diffusion. 

Global simulations of Type I migration of a planet in 
a radiatively inefficient disk wer e performed by som e au- 
thors, using a grid-based code jKlev fc Cridal 120081) and 
a sm oothed particle hydrodynamic code ( Avliffe fc Bate! 
120101 ) . They showed that a planet moves outward in an 
optically thick disk. How ever, the local calculation by 
iMorohoshi fc Tanakal (|2003h showed that the one side torque 
approaches to that in the isothermal disk with an increase in 
the efficiency of cooling by radiation, suggesting the inward 
migration of a planet in a radiatively efficient disk. In this 
study, we examine the total torque acting on a planet by 
radiatively efficient disks with various values of the opacity 
and reveal how the total torque depends on the efficiency 
of the cooling by radiation and the physical structure of a 
disk. 

This paper is organized as follows. In the section 2, we 
briefly describe the basic equations and the initial condi- 
tions. The dissipative term due to radiation is included in 
the energy equation. In the section 3, we show the results 
of two-dimensional hydrodynamic simulations, taking into 
account the effect of the cooling due to radiation. We show 
the dependence of the total torque on the cooling efficiency 
by radiation. The direction and the magnitude of the migra- 
tion of a planet change depending on the thermal structure 
of a disk and the cooling efficiency. Summary is given in the 
section 4. 



equations are the mass conservation, the Euler equations, 
and the energy conservation equation with a dissipative term 
due to radiation by small dust particles in the disk. 

We use a cylindrical coordinate where the star is located 
at the center of the coordinate. The mass of the planet is 
much smaller than that of the star and we neglect the indi- 
rect term. The basic equations read 
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where E is the gas surface density, p and E are the verti- 
cally integrated pressure and the specific energy of the gas, 
respectively; v r and v$ are the radial and tangential veloc- 
ities of the gas; k is the Planck mean opacity, <jsb is the 
Stefan-Boltzmann constant, T is the temperature, Ti n i is 
the unperturbed initial temperature, and <3? is the gravita- 
tional potential. We consider a less massive disk and neglect 
the self-gravity of the gas. The gravitational potential of the 
star and the planet is given by 

GM e GM p 
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where Mq is the mass of the star, M p is the mass of the 
planet, ip is the angle between the r- vector and the r p -vector, 
and Hp is the scale height of the disk at the location of the 
planet. The scale height is given by 
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(6) 



2 BASIC EQUATIONS AND NUMERICAL 
METHOD 

2.1 Basic Equations 

A planet excites density waves in a protoplanetary disk 
and changes the density distribution of a disk. We exam- 
ine the torque acting on a planet by a radiatively efficient 
disk. In our study, a protoplanetary disk is assumed to be 
optically thin because dust particles are expected to de- 
plet e owing to the formation of planetesimals and plan- 
ets |Morohoshi fc Ta naka 200^). We consider a planet with 
5 Earth masses in an optically thin disk. A planet rotates 
around a solar mass star in a fixed circular orbit. The dis- 
tance of the planet from the star is 15 AU and denoted by 
r p . 

The problem is limited to two-dimensional flow, where 
all physical quantities (e.g., the surface density) depend on 
r and 0, where r is the distance from the star and 8 is the 
angle between the x-axis and the position vector. Governing 



where c p and f2 p are, respectively, the isothermal sound ve- 
locity and the Keplerian angular velocity at r p . The smooth- 
ing length parameter, e, is introduced i n order to include the 
effect of the scale height of the disk. iTanaka et al.1 (|2002h 
pointed out that the gravity of disk gas on a planet in a 
two-dimensional disk is greatly enhanced, compared with 
that in a three-dimensional disk. In this study, e is set to be 
0.2 for a planet with 5 Earth masses (See Appendix). 

iMuto fc Inutsukal i|2009T ) made the linear calculation 
and found that the one side torque exerted on a planet by a 
viscous disk converges to a well defined value as long as the 
viscosity is small. We neglect viscosity of the disk gas, even 
though viscosity of the gas prevents the corotation torque 
from saturation. We further assume that the temperature 
of the disk gas is the same as that of the dust particles. 
The dissipation term in the energy equation is derived as- 
suming the local thermal equilibrium in an optically thin 
disk. The efficiency of the energy transport by radiation is 
mainly determined by the Planck mean opacity. When a 
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large number of small dust particles absorb light in a disk, a 
disk becomes opaque. Following the formation of large par- 
ticles by coalescence of the small dust particles, the opacity 
of a disk decreases. The evolution of the opacity strongly 
depends on the growth of small dust particles. Collisions 
of dust particles with large impact energies might result in 
fragmentation of the particles. A large number of small dust 
particles produced by the fragmentation increase the opac- 
ity of the disk. The opacity does not decrease and even in- 
creases when b r ittle d ust particles are considered. However, 
iBirnetiel etafl (|200Sh showed that the opacity significantly 
decreases in a million years if dust particles coalesce when 
the collision velocities of the particles are less than 10 m/s. 

The specific energy of the gas is given with the pressure 
and the kinetic energy: 



7—1 2 x ' 



(7) 



where 7 is t he ratio o f the s pecific heats at constant pressure 
and volume. iLi et al.l (|2000T l obtained 7 of a two-dimensional 
disk as 7 = 4/3. The gas is assumed to be ideal and the 
equation of state is given by 

E/cbT 
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where fee is the Boltzmann constant, mn is the mass of a 
hydrogen atom, and [i is the mean molecular weight of the 
gas. We set up /i =2.34. 

The angular momentum of the disk gas is transferred 
to the planet. The transfer rate of the angular momentum 
from the disk gas at r to the planet is given by 
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By integrating the torque density over the radial distance, 
we obtain the total torque acting on the planet from the 
disk: 



IVdr, 



(10) 



where r max and r m i n are, respectively, the outer and inner 
radii of the disk. 

For the later convenience, we write all the quantities in a 
non-dimensional form using t he unit length r p , the unit mass 
Mq, and the unit time Qp 1 . iTanaka et all (|2002f ) made the 
linear calculation of the gravitational interactions between 
a planet and an isothermal disk. They obtained the torque 
on the planet exerted by the disk and found that the torque 
has the magnitude of r = (M p /M e ) 2 (r p tt p /c p ) 2 T< p r p Q,l, 
where £ p is the surface density at r p . In this study, the 
total torque and the torque density are normalized by To and 
To/Vp, respectively. The normalized quantities are denoted 
with a tilde (e.g., T). 



2.2 Numerical Method 

We consider a disk whose inner and outer radii are given 
by r m i n =6 AU and r max =24 AU. We use an equidistant 
grid in r and 6 directions with th e resolution of 576 x 3072 . 
The killing boundary conditions (|de Val-Borro et al.ll2006T l. 
where all components are relaxed towards their initial state, 
are used in the inner zone (6 AU to 7.5 AU) and in the 



outer zone (21 AU and 24 AU) to avoid wave reflections 
from the boundaries. All the quantities in the inner and 
outer boundaries are always fixed to be the initial values. 

Initially the temperature and the surface density distri- 
butions of the disk are given by 



T ini = T p f- q 
and 

Sini — Spf" P 



(11) 
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respectively, where T p is the temperature at the location of 
the planet and set to be 72 K. The scale height of the disk, 
H, is given by 
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The surface density, E p , is determined such that the total 
mass of the disk becomes the same as th at of the standard 
disk model, 5 x 10~ 3 M© (|Havashilll98lh . The initial rota- 
tional velocity of the gas is calculated from the force balance 
in the r-direction: the gravity of a star, the centrifugal force, 
and the radial pressure gradient, while the initial radial ve- 
locity is zero. Using obtained spectr al energy distributions 
of disks around young stellar objects. iKitamura et al.l (|2002l ) 
showed that the power-law indices of the temperature and 
surface density distributions of disks range from 0.5 to 0.7 
and from to 1.0, respectively, that is, 0.5 < q < 0.7 and 
< p < 1. We calculate the corotation torque and the Lind- 
blad torque acting on the planet by the disk with various 
combinations of the two parameters, p and q, which covers 
the parameter space suggested by the observations. 

We develop a two dimensional global hydrodynamic 
program with the gravitational forces of a star and a planet 
and a dissipative term by radiation. The basic equations are 
solved simultaneously using the finite volume method with 
an operator splitting procedure. The source terms are com- 
puted with a second order Runge-Kutta scheme, while the 
advection terms are calculated with a second ord er MUSCL- 
Hancock scheme and an exact Riemann solver (|Torol [l999l : 
llnaba et alll2005l ). The isothermal Riemann solver is used 
to calculate the total torque on a planet in an isothermal 
disk. It is noted that a disk with the fixed temperature dis- 
tribution is hereafter called an isothermal disk in our study. 

Before calculating the torque on the planet in a disk, let 
us introduce the cooling time, t cao \, given by c\r /16asBT? ni n, 
where cy is the specific heat at constant volume. Further- 
more, we introduce the compression time, t comp . A gas el- 
ement is compressed while it approaches the planet from 
its fron t on the horseshoe orbit. IPaardekooper fc PapaloizovJ 
(|2009bl ) investigated the motion of gas on a horseshoe orbit 
of a low mass planet and obtained the width of the horse- 
shoe region. The half width of the horseshoe region is given 
by 



f s = 2, 



Mp 

Mr, 



(14) 



The velocity of a gas element on the horseshoe orbit is ap- 
proximated by 1.5fi p r s . Assuming that a gas element is com- 
pressed while crossing the horseshoe region (~ 2r s ), we have 
tcomp = 4/ (3fi p ). We define the cooling efficiency by the ra- 
tio of the compression time to the cooling time: 
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C = 



64agBgw 
3cvf2 P 



(15) 



The cooling efficiency indicates how efficient the heat is dis- 
sipated from the gas element while the gas element ap- 
proaches to the planet an d undergoes the compression. 
iMorohoshi fc Tanakal (|2003f ) estimated the disk opacity to 
be 0.0075cm 2 /g when a disk has dust particles with 1cm- 
radius. The cooling efficiency of the disk is given by 



C = 18 



0.0075cm 2 /g 



(16) 



In an adiabatic disk, £ = 0. A disk approaches an isothermal 
disk with an increase in ( due to the efficient cooling of the 
disk. 



3 RESULTS 

3.1 The torque density exerted on a planet 

We study the gravitational interactions between the planet 
and a disk numerically. The gravity of the planet generates 
density waves in a disk inside and outside the orbit of the 
planet. The inner and outer trailing density waves pull the 
planet gravitationally. The positive and negative torques by 
the inner and outer density waves, respectively, act on the 
planet. In a standard disk, the gravity of the outer wave is 
stronger than that of the inner wave because the location 
of the outer density wave is closer to the planet due to the 
pressur e gradient o f a disk, leading to the negative Lindblad 
torque (|Wardlll997l ). The negative Lindblad torque decreases 
the angular momentum of the planet. 

There is another torque exerted on the planet by 
a gas element on t he ho rs eshoe orbit of the planet 
JBaruteau . 



Massetl 120081 ; iPaardekooper fc Papaloizoul 
l2009bh . When a gas element on the horseshoe orbit ap- 
proaches the planet, the angular momentum is exchanged 
between the gas and the planet. This causes the positive 
or negative corotation torque on the planet, depending 
on the thermal structure of the disk. The corotation 
torqu e is dependent on t he en tropy distribution of the 
disk (B aruteau fc Massetl 120081 ). The entropy of a gas 
element in an adiabatic disk is conserved while it moves 
on the horseshoe orbit of a planet. In an adiabatic disk 
with the negative entropy gradient, a gas element on the 
horseshoe orbit is compressed when its orbit is changed 
from a slower orbit to a faster orbit, increasing the density 
of the gas element. The positive torque is exerted on the 
planet by the gas element. Therefore, the corotation torque 
is dependent on the entropy distribution of the disk. The 
initial entropy distribution is given by S — p/E 7 oc r A , 
where A is (7 — l)p — q. 

It was suggested that the corotation torque greatly 
decreases after the synodic period of the planet be- 
cause the entropy of the gas in the horseshoe region of 
the planet tends to become uniform (jBaruteau fc Massetl 
20081 ; IPaardekooper fc Papaloizoul 12008 , 12009 j ; iKlev et all 
20091 '). This is called the saturation o f the coro- 
tation torque. Paa rdckoopcr & Pap aloizoul (|2008l ) and 
IPaardekooper fc Papaloizoul ( 2009al ) showed that the coro- 



tation torque does not saturate in viscous disks because fresh 



gas is supplied into the horseshoe region from the outer re- 
gion of the disk by viscosity. We also find the saturation of 
the corotation torque after the synodic period of the planet 
(about 30 rotations of the planet) because viscosity is not 
included in this study. If we include viscosity, we might not 
find the saturation of the corotation torque. Since we would 
like to study the effect of energy dissipation by radiation on 
the corotation torque in this study, we stop the simulations 
before the corotation torque saturates. Disks reach a tenta- 
tive steady state in about ten rotations through successive 
gravitational interactions with a planet. 

We first consider the planet in the isothermal and adi- 
abatic disks to find the difference of the torques in the ide- 
alized disks. The disks have the initial density and temper- 
ature distributions with p — 0.8 and q = 0.7, which gives 
the negative entropy gradient, A = —0.43. Fig. [1] shows the 
radial distribution of the torque density, the equation ([9]), 
exerted on the planet by the isothermal and adiabatic disks 
at t/tp = 20, where t p is the rotational period of the planet. 
The solid and dashed curves correspond to the torque den- 
sities exerted on the planet by the adiabatic and isother- 
mal disks, respectively. The large corotation torque is found 
in the adiabatic disk. The corotation point is located at 
f c = 1 — 6 x 10 -3 , where f c is the non-dimensional distance 
of the corotation point from the star. The torque density in 
the adiabatic disk becomes maximum at f ~ 1 — 5 x 10 -3 . 
The corotation torque appears in the adiabatic disk due to 
the enhancement of the gas density on the horseshoe or- 
bit. Furthermore, the torque density has a local maximum 
and a local minimum at f ~ 1 — 6 x 10 -2 and 1 + 6 x 10 -2 , 
respectively. The distances from the planet to the local max- 
imum and minimum are approximated by the scale height 
of the disk. The local maximum and minimum of the torque 
density correspond to the Lindblad torque. In the isother- 
mal disk, the inner one side torque is 8.6, while the outer 
one side torque is —11.1. As a result, the total torque ex- 
erted on the planet by the isothermal disk becomes negative 
(r = —2.5). The planet migrates inward toward the star in 
the isothermal disk. On the other hand, the large corotation 
torque found in the adiabatic disk becomes larger than the 
negative Lindblad torque and, as a result, the total torque 
becomes positive (F = 0.3). The planet in the adiabatic disk 
migrates outward. 

We further perform a simulation and calculate the 
torque density exerted on the planet by the disk that dissi- 
pates the energy by radiation to examine the effect of the en- 
ergy dissipation on the planetary migration process. We as- 
sume that the disk contains a number of 1 cm dust particles 
responsible for the energy dissipation. We adopt £ = 18 as in 
the equation (|16[) . Hereafter we call a disk with the energy 
dissipation by radiation a radiative disk. Fig.[2]shows the ra- 
dial distributions of the torque density exerted on the planet 
by the adiabatic disk and the radiative disk at t/tp = 20. 
The initial disk structure is the same as that in Fig. [T] The 
solid and dashed curves correspond to the torque densities in 
the adiabatic and radiative disks, respectively. The corota- 
tion torque in the radiative disk becomes smaller than that 
of the adiabatic disk because the adiabatic compression gets 
weaken by the energy dissipation. As shown in Fig. [2] the 
magnitude of the Lindblad torque density in the radiative 
disk is larger than that of the adiabatic disk. The torque 
density distribution of the radiative disk becomes similar to 
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that of the isothermal disk. The total torque becomes neg- 
ative (f 1 = —2.5), leading to the inward migration of the 

planet. 

iMorohoshi fc Tanakal (|2003h performed the local shear- 
ing box simulations of the gravitational interactions between 
a planet and a radiative disk. They plotted the contour of 
the gas density near a planet and found the elliptical shape 
of the contour. The direction of the semimajor axis of the 
contour is inclined with respect to the direction toward a 
central star. This results in the larger torque density in the 
vicinity of a planet because the difference in the density dis- 
tribution between the front and rear of the planet is larger. 
The contour of the density distribution in the vicinity of the 
planet is shown in Fig. [3] at t/t p = 20. The panels (a) and 
(b) show the density contours in the adiabatic and radiative 
disks, respectively . Our numerical resul t s are qualitatively 
similar to that of IMorohoshi fc Tanakal (|2003T ). As seen in 
Fig- E] (a), the gas density is enhanced in the inner horseshoe 
orbit region, yielding the larger positive corotation torque in 
the adiabatic disk. The gas density around the planet in the 
radiative disk becomes larger than that in the adiabatic disk. 
Heat generated by the compression is dissipated by radia- 
tion and the gas density is required to increase to sustain the 
hydrostatic equilibrium in the vicinity of the planet. Hence, 
the larger Lindblad torques in the radiative disk than those 
in the adiabatic disk are found in Fig. [2] 

3.2 Dependence of the total torque on the 
radiative efficiency and the planet mass 

We carry out a number of numerical simulations with vari- 
ous values of the radiative efficiency, (. Fig. [4] shows the time 
evolution of the total torque exerted on the planet by the 
disk with p = 0.8 and q = 0.7. The solid, dotted, dashed, 
long-dashed, and dot-dashed curves correspond to the total 
torques in the cases of £ =0, 0.18, 0.9, 1.8, and 18, respec- 
tively. The total torque increases with time at the begin- 
ning of the simulations and reaches a steady state around 
t/t p = 10. Note that the saturation of the corotation torque 
occurs after the synodic period (t/t p ~ 30) if we continue the 
simulation. The total torque decreases with an increase in 
the radiative efficiency. The corotation torque becomes small 
when the heat generated by the compression is dissipated 
by radiation effectively. The contribution of the Lindblad 
torque becomes dominant and the sign of the total torque is 
changed from positive to negative with the increasing radia- 
tive efficiency. Except in the adiabatic disk, the total torque 
becomes negative and the planet moves inward toward the 
central star in the radiative disks. In Fig. O the total torque 
is shown as a function of £ for the disk with p = 0.8 and 
q — 0.7. The total torque decreases with increasing ( and 
approaches that by the isothermal disk. The total torque 
acting on the planet by the radiative disk with f > 9 is 
nearly the same with that by the isothermal disk. 

Fig. [6] shows the time evolutions of the total torques 
exerted on the planet by the radiative disks with some ra- 
diative efficiencies, In the panel (a), we use the disk with 
p — 0.5 and q = 1.0. The solid, dotted, dashed, and dot- 
dashed curves correspond to the total torques exerted on the 
planet by the disks with ( — 0, 0.18, 1.8, and 18, respectively. 
The disk has the negative entropy gradient and the power- 
law index of the entropy distribution is given by A = —0.83. 



This is smaller than that of the disk used in Fig. [4j exerting 
the larger positive corotation torque on the planet. The to- 
tal torques by the disks with £ = (the adiabatic disk) and 
( — 0.18 become positive due to the large positive corotation 
torque. The corotation torque decreases with the increasing 
radiative efficiency £. The disk with £ > 1.8 exerts the neg- 
ative torque on the planet. In the panel (b), we use the disk 
with p = 1.0 and q — 0.5. The solid, dotted, dashed, and 
dot-dashed curves correspond to the total torques by the 
disks with £ = 0, 0.18, 0.9, and 1.8, respectively. The power- 
law index of the entropy distribution is A = —0.16, which is 
larger than that of the disk in Fig. 2] The gas on the horse- 
shoe orbit of the planet is compressed weakly in the disk 
with large A, making the corotation torque small. The total 
torque becomes negative even in the adiabatic disk. 

The total torque is dependent on the radiative efficiency, 
(, and the power-law index of the entropy distribution, A, 
as shown in Fig. [7| We consider the disks with p + q = 1.0 
and 1.5. The total torques at t/t p = 20 are plotted. The 
filled circles, the open squares, the filled diamonds, and the 
open triangles correspond to the total torques on the planet 
exerted by the disks with the radiative efficiencies of £ = 0, 
0.18, 1.8, and 18, respectively. The fluctuations of the total 
torques are due to the different choices of p and q. Table [1] 
shows the list of the parameter values of p and q a s wel l 
as the thermal state of the disk. iBaruteau fc Massed (|2008T l 
showed that the corotation torque increases with decreasing 
p+q as well as decreasing A. We confirm their results and find 
that the total torque exerted by the disks with p + q — 1.0 is 
larger than that by the disks with p + q = 1.5. The Lindblad 
torque is always negative in disks with the negative pressure 
gradient. On the other hand, the corotation torque increases 
to be positive with decreasing A. The total torque becomes 
positive in the adiabatic disks with A < —0.4. The planet 
migrates outward in the adiabatic disks with A < —0.4. 

The total torque decreases when the energy dissipation 
is included. The total torque on the planet by the disk with 
( — 0.18 and A = —0.83 becomes the half of that by the adi- 
abatic disk. The energy dissipation by radiation decreases 
the corotation torque, resulting in the decrease in the to- 
tal torque. In the disks with £ = 0.18, the sign of the total 
torque changes at A ~ —0.5. When ( > 1.8, the total torque 
is always negative. The planet migrates inward toward the 
central star even in the disk, of which the amount of dust 
particles is one-tenth of that adopted in the standard disk. 
The total torque approaches that exerted by the isother- 
mal disk with increasing (. The total torques exerted on the 
planet by the disks with ( — 1.8 and 18 are weakly depen- 
dent on the power-law index of the entropy distribution. The 
contribution of the corotation torque to the total torque be- 
comes small. Even if the total torque by the adiabatic disk 
becomes positive due to the large positive corotation torque 
by the adiabatic compression of gas on the horseshoe orbit, 
it becomes negative when we take into account the energy 
dissipation in the disk. The planet in the disk with very 
small A and ( might migrate outward. 

iKitamura et all (|2002T ) observed the thermal emission of 
dust in protoplanetary disks around T tauri stars in Taurus. 
They derived the disk properties of the surface density and 
temperature distributions. The power-law indices p and q 
range from to 1 and from 0.5 to 0.7, respectively. If we 
only consider disks with the age over 10 6 yr, the power-law 
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indices occupy the narrower domains of < p < 0.6 and 
0.5 < q < 0.6. Then the range of the power-law index A 
of the entropy distribution becomes —0.6 < A < —0.3. The 
disks with —0.6 < A < —0.3 exert the positive torque on 
the planet when the dust opacity of the disk is smaller than 
1/100 of the standard dust opacity (see the equation HU)). 
The depletion of dust particles results in the outward or very 
slow inward migration of the planet in the observed disks. 

Once we obtain the total torque, F, exerted on the 
planet by the disk, the r adial migration vel ocity of the 
planet, r p , is calculated as (|Tanaka et alj|2002f ) 

(t)[*W <"» 

where L p is the angular momentum of the planet given by 
L p = Mp(GM r p ) 1/2 and M® is the Earth mass. The neg- 
ative and positive radial velocities of r p correspond to the 
inward and outward migration of the planet, respectively. As 
shown in Fig. [Jj the total torque is on the order of —1 in the 
radiative disk with the large £ and the planet moves about 
1 AU in 10 yr. The large migration velocity prohibits the 
planet from growing. As shown in Fig. [Jj the total torques 
in some radiative disks with small £ happen to be on the 
order of 0.1. The planet embedded in these disks migrates 
slowly and might be able to grow further. 

Finally we study the dependence of the total torque on 
the planet mass. The total torques exerted on the planet 
with 20 Earth masses by the disks with various A and 
£ at t/tp — 10 are shown in Fig. [8] The total torque 
quickly reaches a steady state because the width of the 
horseshoe orbit of the planet is wide and it takes short 
time for a gas element on the horseshoe orbit to approach 
the planet. We find the width of the horseshoe orbit of 
the planet with 20 Earth masses is tw o times as wide as 
that of the planet with 5 Earth ma sses l|Masset et al.ll2006l ; 
iPaardekooper fc PapaIoizo"u]|2009bT l. We use the disks with 
p + q = 1.0 and 1.5. The filled circle, the open square, 
the filled diamond, and the open triangle correspond to 
the total torques by the disks with ( = 0, 0.18, 1.8, and 
18, respectively. The sign of the total torque on the planet 
with 20 Earth masses in the adiabatic disk changes around 
A = —0.4, which coincides with the case of the planet with 
5 Earth masses as shown in Fig. [7] The zero total torques 
on the planets with 5 and 20 Earth masses are found at the 
same value of A ~ — 0.6 in the disks with £ = 0.18 as well. 
iKlev fc Cridal (|2008T l studied Type I migration of the planet 
with 20 Earth masses in the disk with p — 0.5 and q = 1.6. 
The disk model has the very steep entropy gradient with 
A = —1.4. Due to the steep entropy gradient, the planet 
moves outward in the disk even with large opacity. 

The linear analysis of the gravitational interactions be- 
tween a planet and an isothermal disk showed that the to- 
tal torque is propor tional to the square of the planet mass 
|Tanaka et al.ll2002i ). Fig. [5] shows the total torques for var- 
ious planet masses in the disk with p = 0.5 and q — 1.0. 
Although the linear analysis suggested the constant magni- 
tude, one can see from Fig. [5] that it decreases with the 
increasing planet mass. This discrepancy would orig inate 
from the non- linear effect. iMorohoshi fc Tanakal (|2003l ) also 
showed the non-linear effect when large planets are consid- 
ered. As compared between Figs. [Jj and [8] the total torque 
exerted on the planet with 20 Earth masses becomes some- 



what smaller for the positive values and larger for the neg- 
ative values than that on the planet with 5 Earth masses in 
whole. 



4 SUMMARY 

A planet exerts the gravity on gas in a disk and changes 
the density distribution of the gas. Trailing density waves 
are generated inside and outside the orbit of the planet. 
The outer density wave pulls the planet gravitationally and 
exerts the negative torque on the planet, while the posi- 
tive torque is exerted on the planet by the inner density 
wave. The magnitude of the gravitational force by the outer 
wave on the planet is usually larger than that by the inner 
wave due to the negative pressure gradient. The sum of the 
torques exerted by the outer and inner waves become s neg- 
ative , leading to the inward migration of the planet (|Wardl 
Il997t ). The torque exerted on the planet by the density waves 
is called the Lindblad torque. A gas element on the horseshoe 
orbit of a planet exerts torque on the planet as well. The an- 
gular momentum is exchanged between the planet and the 
gas element on the horseshoe orbit while the gas element 
moves from a slow orbit to a fast orbit. The gas element is 
compressed when the entropy gradient of the disk is nega- 
tive. The compressed gas element increases the density and 
exerts the stronger gravity on the planet. The torque exerted 
on the planet by the gas on the horseshoe orbit is called the 
corotation torque. The corotation torque becomes positive, 
leading to the outward migration of the planet, when the 
entropy gradient of the disk is negative. The energy dissipa- 
tion has an important role to determine the magnitude of the 
corotation torque. The density of a gas element on the horse- 
shoe orbit increases the most in an adiabatic disk, in which 
a gas element conserves the entropy during the change of 
the orbit. The energy dissipation reduces the entropy of the 
gas element and makes the increment of the density small. 
The corotation torque decreases when the energy dissipation 
is included. We include the energy dissipation by radiation 
in this study. The sum of the Lindblad and the corotation 
torques, the total torque, determines the direction and the 
magnitude of the migration of the planet. 

The total torques exerted on the planet by the isother- 
mal and adiabatic disks were first examined. We find that 
the total torque becomes negative in the isothermal disk. 
The planet in the isothermal disk migrates inward toward 
the central star. On the other hand, the total torque becomes 
positive in the adiabatic disk when the entropy distribution 
of the disk has a steep negative slope. Gas elements on the 
horseshoe orbit are compressed adiabatically. The increased 
density during the change of the orbit from a slow orbit to 
a fast orbit enhances the corotation torque. The corotation 
torque saturates in a synodic period of the planet because 
viscosity is not included in the present study. The magni- 
tude of the corotation torque is dependent on the value of 
the power-law index of the entropy distribution, A. The gas 
disks with the negative A (the positive A) exert the posi- 
tive (negative) corotation torques on the planet. The total 
torque exerted on the planet by the adiabatic disk becomes 
positive when A < —0.4 and the planet migrates outward in 
the disk. 

We have further made a number of numerical Simula- 
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tions of the gravitational interactions between the planet 
and the disk, taking into account the energy dissipation by 
radiation. The density distribution of gas around the planet 
in the radiative disk is different from that in the adiabatic 
disk. More gas pressure is required to balance the gravity 
of the planet in the radiative disk because of the energy 
dissipation. The density around the planet in the radiative 
disk is larger than that in the adiabatic disk. The elliptical 
shape of an isoline of the density, whose semimajor axis is in- 
clined with respect to the direction toward the central star, 
is formed near the planet. It is unclear why the isolines are 
more inclined in the radiative disk. Because of the inclined 
isoline of the density, the larger Lind blad torque is exerted 
on th e planet in the radiative disk (|Morohoshi fc Tanakal 
120031 ). On the other hand, the corotation torque decreases 
with an increase in the energy dissipation. We consider the 
energy dissipation by radiation. The efficiency of the energy 
dissipation, (, is defined by the ratio of the two time periods: 
the time period for a gas element on the horseshoe orbit to 
change the orbit and the time period for a gas element to 
lose the energy by radiation. Energy dissipates effectively in 
the disk with large The corotation torque decreases with 
an increase in leading to the negative total torque. The 
sign of the total torque is dependent on The total torque 
invariably becomes negative in an optically thin disk with 
the reasonable range of A. Even though the sign of the total 
torque is negative, the magnitude of the total torque might 
become small when the disk has small (,. 

The total torque is also dependent on the mass of the 
planet. We examine the total torques exerted on the planets 
with 5 and 20 Earth masses. The total torque is approxi- 
mately proportional to the square of the mass of the planet 
as shown by the linear analysis. With an increase in A, the 
total torques exerted on the planets with 5 and 20 Earth 
masses by the adiabatic disk decrease and change the sign 
at A ~ —0.4. The direction of Type I migration is little 
sensitive to the planet mass. 
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APPENDIX A: DERIVING THE SOFTENING 
PARAMETER 

A planet interacts with a gas disk gravitationally. We con- 
sider a two-dimensional gas disk, in which the gas exists only 
in a plane. The planet in the same plane induces the larger 
gravity on the gas owing to the close distance. We introduce 
a softening parameter to incorporate the effect of the scale 
height of the disk. The softening parameter is used to reduce 
the gravity of the planet on the gas in the disk. 

We assume the hydrodynamic equilibrium of the gas in 
the vertical direction. The velocity of the gas is independent 
of z. The temperature of the gas is assumed to be constant 
in the ^-direction because an optically thin disk is consid- 
ered in this study. The gravity of the planet on the three- 



8 K. Yamada and S. Inaba 



200 
100 

r r 

-100 



0.8 0.9 1^ 1.1 1.2 

~r 

Figure 1. The torque density acting on the planet with 5 Earth masses at t/t p = 20, where t p is the rotational period of the planet. The 
disk has the density and temperature distributions with p = 0.8 and q = 0.7, which gives the negative power-law index (A = —0.43) of the 
entropy distribution. The torque densities by the adiabatic and isothermal disks are shown by the solid and dashed curves, respectively. 
The positive corotation torque appears due to the negative gradient of the entropy distribution in the adiabatic disk. The torque density 
becomes zero in a far place from the planet due to the weak gravitational interactions. We plot the torque density in the range of 
0.8 < f < 1.2, even though we calculate the gravitational interactions between the planet and the gas disk with 0.4 < f < 1.6. 
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Figure 2. The same as described in the legend to Fig.[l] but the dashed curve is the torque density by the radiative disk with £ = 18. 



dimensional gas is expressed by pV$3D, where $3d and the 
density p, respectively, are given by 



$3D = 



and 



P = 



exp 



(Al) 



(A2) 



We integrate the gravity of the planet over z as 

pV$ 3 od2. (A3) 



Note that V does not include the partial derivative with 
respect to z and d/dz. 

The softening parameter, e, is defined implicitly by the 
following equation: 



EV- 



GM 



pV<I>3 D d2. 



(A4) 



v/|r-r p p + e^ 

Assuming the constant scale height, H p , near the orbit of 
the planet, we have the following relation: 

1 = J_ r°expH,V). 
y/y 2 + e 2 V* Jo 

where y and x are defined by |r — r p \/H p and z/\r — r p | 
respectively. The softening parameter, e, is calculated as 



-da;, 



(A5) 



exp(— y 2 x 2 ) 
VI + x 2 



Ax 



(A6) 



The softening parameter is shown in Fig. IA1I The soft- 
ening parameter gradually increases with an increase in y. 
In our study, e is set to be 0.2 and 0.3 for a planet with 
5 Earth masses and 20 Earth masses, respectively. This soft- 
ening parameter roughly corresponds to the Bondi radius 
(y — 1.6 x 10~ 2 ) of the planet with 5 Earth masses and the 
Bondi radius (y = 7.0 x 10 -2 ) of the planet with 20 Earth 
masses. 
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Figure 3. The contour plot of the gas surface density around the planet at t/t p = 20 in (a) the adiabatic disk and (b) the radiative 
disk with f = 18. The initial surface density is subtracted from the surface density for the better contrast. The densities are normalized 
by the initial density at r p . The planet is located at f = 1 and 9 = ir. The gas density near the planet is enhanced in the radiative disk, 
compared with that in the adiabatic disk because heat is quickly radiated and the gas density is required to increase to maintain the 
hydrostatic equilibrium. 
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Figure 4. The time evolution of the total torque exerted on the planet by the disks with the various radiative efficiencies: f = (solid), 
0.18 (dotted), 0.9 (dashed), 1.8 (long-dashed), and 18 (dot-dashed). The disk with f = corresponds to the adiabatic disk. The total 
torque becomes positive in the adiabatic disk. With an increase in f , the disk loses the energy effectively and approaches the isothermal 
disk. Initially all the disks have the power-law distribution of the density and the temperature with p = 0.8 and q = 0.7, respectively. 
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Figure 5. The total torque exerted on the planet by the disk with the various cooling efficiencies at t/t p = 20. The dotted line shows 
the total torque exerted by the isothermal disk, which has the fixed temperature distribution with the power-law index of q = 0.7. The 
total torques exerted by the disks with £ > 9 is nearly the same as that by the isothermal disk. 




Figure 6. The time evolutions of the total torques exerted on the planet by the two different disks with (a) p = 0.5 and q = 1.0 and (b) 
p = 1.0 and q = 0.5. The solid, dotted, dashed, and dot-dashed lines correspond to f = 0, 0.18, 1.8, and 18 in the panel (a) and £ = 0, 
0.18, 0.9, and 1.8 in the panel (b), respectively. In the panel (a), the disks with £ = and 0.18 exert the positive torque on the planet. 
In the panel (b), the total torque becomes negative even in the adiabatic disk. 
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Figure 7. The total torque exerted on the planet with 5 Earth masses by the disk as a function of the power-law index of the entropy, 
A, in the cases of £ = (filled circle), 0.18 (open square), 1.8 (filled diamond), and 18 (open triangle). 




Figure 8. The same as described in the legend to Fig. [7] but the planet with 20 Earth masses. 
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Figure 9. The total torque as a function of the planet mass embedded in the disk with p = 0.5 and q = 1.0. The filled circle and 
diamond correspond to the cases of C, = and 18, respectively. The linear analysis showed that the total torque is proportional to the 
square of the planet mass. The total torque divided by the square of the planet mass decreases with an increase in the planet mass. This 
is the non-linear effect. 
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Table 1. The parameter values and the thermal state of the disk used for the simulations 



Model number 
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adiabatic 
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adiabatic 
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adiabatic 
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adiabatic 
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-0.3 


adiabatic and radiative 
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adiabatic and radiative 
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adiabatic and radiative 
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adiabatic and radiative 
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Figure Al. The softening parameter as a function of y. 



